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New nuclear structure function data from Jefferson Lab covering the higher x and lower Q 2 
regime make it possible to extract the higher order Fi moments for iron and deuterium at low 
four-momentum transfer squared Q 2 . These moments allow for an experimental investigation of 
the nuclear momentum sum rule and a direct comparison of the non-singlet nucleon moment with 
Lattice QCD results. 
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I. INTRODUCTION 

Nuclear effects in lepton-nucleus scattering have been 
extensively studied, both experimentally and theoreti- 
cally, over the last few decades. For recent reviews, see 
Refs. PIQ- The body of available data provides clearcut 
evidence that the nucleus can not be simply described 
as a collection of nucleons on mass shell. For example, 
the study of nuclear structure functions led to the dis- 
covery of the "EMC effect" where it was found that the 
quark distribution inside the nucleus differs from that 
of a free nucleon. The availability of experimental in- 
formation on the Q 2 dependence of the moments of the 
nuclear structure function F^(x, Q 2 ) has stimulated the- 
oretical analyses of meson exchange contributions and 
off-shell effects in nuclei, sometimes showing sizeable 
deviations from predictions of simple convolution mod- 
els 0, 0, IU El- In this discussion, A is the mass num- 
ber, Q 2 is the four-momentum transfer squared in the 
lcpton-nucleon inclusive scattering process, and x is the 
Bjorken scaling variable, with < x < 1 for the proton, 
< x < Ma/M p ps A for a nucleus. 

Previous nuclear structure function moment analyses 
have relied on moment data extracted from several ex- 
periments carried out at CERN 0, and SLAC 0, E3 
using 56 Fe and 2 H targets. The experimental values of 
Cornwall-Norton moments, M n (A, Q 2 ), require precision 
measurements of structure functions covering large inter- 
vals of x, Q 2 , and A, since: 

M n (A,Q 2 )= / dxFf(x,Q 2 ) x n - 2 . (1) 
Jo 

Here, n is an integer defining the order of the moments. 
We note that the n = 2 moment can be related to the 
familiar momentum sum rule, which must be less than 
unity for the nucleon. Asymptotically, QCD predicts the 
fraction of the nucleon momentum carried by the quarks 
to be (1 + 16/3/) -1 , where / is the number of quark 
flavors [Tl| . 

Until recently, the set of experimental data at large 



x was rather poor, and thus the evaluation of the mo- 
ments was correspondingly imprecise, especially for large 
n. Typically, data were obtained in the deep inelastic 
scattering regime at moderate to small values of x and 
larger values of Q 2 . One can see immediately from Eq.^ 
that, as n increases, larger x data will increasingly dom- 
inate the moments. Additionally, at lower values of Q 2 , 
the structure function is larger in the higher x region 
and dominates even the lower order moments. Moreover, 
nuclear structure effects are expected to show up most 
clearly at large values of x 0. 

Recently, data have become available from new exper- 
iments at Jefferson Lab which cover higher x and lower 
Q 2 0,0,0,0, complementing the previous data set. 
These new data make it possible to accurately extract 
the moderate and lower Q 2 moments, and moments to 
higher orders. We report here results from a new extrac- 
tion of the F% structure function moments for iron and 
deuterium and compare to proton data. 



II. EXPERIMENT 

Sample spectra used for the extraction of the moments 
are shown in Fig. ^ for deuterium at Q 2 — 4.5 and iron 
at 5 GeV 2 . As noted above, the calculation of the mo- 
ment of a structure function requires data covering the 
whole range in x from to « A at a fixed Q 2 . The struc- 
ture function data used in this analysis were obtained 
in experiments at SLACjllJl|L CERN HEH, Fermi- 
lab [M E3 and JLab [IlIllEl El The Q 2 values 
where the best coverage in x was available were selected. 
In some cases, the data were obtained not at exactly the 
same Q 2 value. In these cases, a small range in Q 2 , vary- 
ing from 0.01 GeV 2 at low Q 2 to 0.5 GeV 2 at high Q 2 
was utilized. The variations of the structure function 
over such ranges were smaller than 2%. 

As illustrated in Fig. ^ the data sets still do not 
cover the full range in x, some extrapolations were nec- 
essary. Between data sets, two methods were utilized, 
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FIG. 1: (Color online) Example structure function data for 
deuterium (top) and iron (bottom) at Q 2 = 4.5 GeV 2 and 5 
GeV 2 , respectively. 



a spline fit and a simple linear extrapolation. Moments 
obtained in such cases agreed within 2%. To extrapolate 
to x — > a parameterization from NMC j2^ was used 
for Q 2 > 2 GeV 2 and a linear extrapolation was used for 
lower Q 2 data. The uncertainty introduced by this ex- 
trapolation becomes negligible for higher order moments. 
The extrapolation to x — > A, while negligible for n = 2, 
becomes important for the higher moments. The data 
used in this region were obtained at SLAC and JLab and 
the coverage in x is sufficient for most Q 2 and n values. 
The uncertainty in the moments due to the extrapolation 
to x = A is less than 1% for n = 2, around 3% for n = 4, 
up to 6% for n = 6, and up to 20% for n = 8. The high- 
est x quasielastic and elastic contributions , im portant for 
low Q 2 , were calculated according to [23|, |24j and added 
to the moments. 

In the near future, the extrapolations to x — ► and 
x — > A can be improved with new data coming from Jef- 
ferson Lab experiments [23 . [2a . 1271 [28l t29l l30l t3l| . These 
experiments have already acquired data and results will 
become available over the next few years. These newer 
data will allow for moments to be obtained over an ex- 
panded range in Q 2 , and for several additional nuclei, 
including 3 He and 4 He. 



III. RESULTS 

Tables [|] and [D] show the Cornwall-Norton moments 
for deuterium and iron. The uncertainties include pub- 



Q 2 

(GeV 2 ) 


n = 2 


n = 4 


n — 6 


n = 8 


0.05 


.481±.481 


.807±.400 


2.3618L.2362 


8.5266L.8527 


0.10 


.407±.204 


.479±.120 


1.0533±.0105 


3.3723L.3372 


0.20 


.320±.080 


.284±.034 


0.3946L.0395 


0.7653L.0765 


0.45 


.296±.021 


.193±.019 


0.2163±.0216 


0.2968L.0359 


0.80 


.220±.011 


.092±.005 


0.0844±.0060 


0.0961L.0103 


1.50 


.180±.009 


.040±.003 


0.0261±.0020 


0.0235L.0033 


2.40 


.169±.008 


.028±.001 


0.0165±.0010 


0.0156L.0008 


3.20 


.162±.008 


.021±.001 


0.0091±.0005 


0.0065L.0003 


4.50 


.165±.008 


.016±.001 


0.0056L.0003 


0.0039L.0002 


5.00 


.161±.008 


.017±.001 


0.0052±.0003 


0.0030±.0002 


7.00 


.163±.008 


.016±.001 


0.0038L.0002 


0.0015L.0001 



TABLE I: Moments of the F2 structure function per nucleon 
for the deuteron. 



lished experimental uncertainties on the structure func- 
tions, the uncertainties due to the finite Q 2 range of 
the data and interpolation procedures, extrapolations to 
low and high x, and the uncertainties in estimating nu- 
clear elastic and quasielastic contributions. The com- 
bined uncertainties are typically 5%, except for low Q 2 
values where the uncertainty in the quasielastic become 
very large, especially for n = 2. At low Q 2 , the higher 
moments become increasingly dominated by the nuclear 
elastic contribution, which is known to better than 5%. 
For the iron n = 6 and n = 8 moments, the intermedi- 
ate Q 2 values have large contributions from the poorly 
known quasielastic contributions at extremely large x val- 
ues, and so these moments are not included. 

There are indications that two-photon exchange cor- 
rections to the electron-nucleon elastic cross section 
might impact the extracted moments |32j . These correc- 
tions appear to be ^6% for elastic e-p scattering (^3% 
for e-n |33j), peaking at large scattering angles. For the 
data included in this analysis, we expect that two-photon 
exchange will contribute at most 2% to the moments, typ- 
ically much less. This is small compared to the experi- 
mental uncertainties, and these effects should partially 
cancel when comparing different nuclei. 

The lower n moments display a very shallow to negli- 
gible Q 2 dependence. In the Operator Product Expan- 
sion (OPE), higher twist effects (interactions between the 
struck quark and other quarks in the electron-nucleon 
scattering process) are expected to manifest a l/Q 2 de- 
pendence in the moment. This is not observed in the 
data, which is somewhat surprising at these low Q 2 values 
where such effects could be large. The asymptotic behav- 
ior of the second moment is ultimately governed by the 
energy-momentum tensor in the OPE and, thus, has no 
Q 2 dependence, as in the quark-parton model ll]. Even 
at the low Q 2 values studied here, the moments display 
this quark-parton model behavior over most of the Q 2 
range. The lower Q 2 moments are dominated by high x 
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(2.0±.l)xl0 9 
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0.40 
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1.00 


.209±.010 


.095±.005 
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— 


1.90 


.166±.008 


.034±.002 


0.0270±0.0015 


.0447±.0058 


2.90 


.174±.009 


.018±.001 


0.0114±0.0010 


.0146±.0063 


5.00 


.158±.008 


.015±.001 


0.0050±0.0004 


.0032±.0006 


6.00 


.164±.008 


.016±.001 


0.0038±0.0002 


.0020±.0004 



TABLE II: Moments of the F2 structure function per nucleon 
for iron. 



resonance regime. Hence, this observation is yet another 
striking manifestation of quark-hadron duality |34| . 

The higher n moments, on the other hand, do display 
an increased Q 2 dependence. These data may therefore 
be used for precision higher twist extractions. However, 
the higher n moments are increasingly dominated by the 
high x, including the elastic and quasi-elastic regimes, 
where the x and Q 2 dependences are less well understood 
in terms of the OPE. 

If nuclear effects are small, the moments for iron can 
also be constructed by adding the proton and neutron 
contributions, extracted from proton and deuteron 
data. To investigate how well this simplified approach 
works, the following simple formula was employed: 

M n (Fe) = Zx M n (p) + (A - Z) x M n {n), (2) 

where M n (n) is taken to be M n (d) — M n (p). Here, 
M n (p), M n (n), and M n (d) refer to the nth moment of the 
proton, neutron, and deuteron, respectively, and Z is the 
atomic number of iron. This is equivalent to extracting 
the iron data as 28 deuterons with a small neutron ex- 
cess contribution. Simple Fermi motion should not yield 
a significant nuclear dependence in the M 2 moment, and 
off-shell effects have been studied 0, @ and are also ex- 
pected to be small for the lowest moment, and on the 
order of 10% for moments up to n = 5 

This procedure is illustrated in Fig. |21 for the sec- 
ond moment, M%. The iron data are shown as squares, 
deuteron data as full circles, proton data as stars. The 
red dashed lines describing the deuteron and the proton 
moments are simple fits to the proton and deuteron data, 
which are then used to calculate the neutron moment, 
Ma(d)— M2(p), and the iron moment as 26 protons and 30 
neutrons, as described in Eq.|2 No additional correction 
was made for nuclear effects or the non-isoscalarity of the 
target. The neutron and iron moments thus calculated 
are shown as blue dashed lines in the figure, while the 
hollow circles show the neutron moments taken directly 
from the difference of deuteron and proton moments. For 
Q 2 > 4 GeV 2 , the ratio of M 2 (Fe)/M 2 (D), normalized 
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FIG. 2: (Color online) The second moment of F2 for proton 
(stars), deuteron (full circles), and iron (squares). The hollow 
circles are the neutron moments taken from the difference of 
deuteron and proton. The red dashed lines are fits to the 
deuteron and proton moments, and the blue dashed lines are 
the neutron and iron moments extracted from these fits using 
the procedure described in the text. 



to the number of nucleons, is 0.99 ±0.05, consistent with 
the value 0.96 j^ij, from a calculation based on Ref. 

The ratios of the measured moments for iron com- 
pared to the moments taken from the deuteron and pro- 
ton moments are shown in Fig. [3] It can be seen that 
these two methods yield the same results within the un- 
certainty. Combining all of the values yields a devia- 
tion of (0.9 ± 2.2)%, or (0.5 ± 2.9)% if we consider only 
Q 2 > 2.5 GeV 2 . This result contradicts interpretations 
of the EMC effect that predict significant modification 
to the the total quark momentum distribution in nu- 
clei. However, it is consistent with other interpretations 
where the total quark momentum is conserved @, H3- 
Here, the data indicate that the integrated iron nucleus 
can be described well as simply being composed of free 
deuterons, with a minimal correction for neutron excess 
in 26p + 30n. It seems the EMC effect is a redistribution 
of quark momentum without any additional momentum 
added by the nuclear environment outside of whatever is 
already present in the deuteron. 

One can also connect the nuclear dependence of the 
quark distributions to the coordinate space parton dis- 
tributions [3^, |3j| . The A dependence of the n = 2 mo- 
ment is then related to the A dependence of the light 
cone distributions at short distance. The fact that the 
data indicate extremely small nuclear effects is consis- 
tent with the result that the A dependence for distances 
less than the inter-nucleon spacing is surprisingly small 
(<2%), due to cancellation between the shadowing, anti- 
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FIG. 3: (Color online) The ratio of the QCD moments for 
iron calculated using iron data to the moments constructed 
using deuterium and proton data shown as a function of Q 2 . 

shadowing, and EMC regions |3Sj . 

We note, further, that the redistribution can be quite 
large, locally. In the structure functions at fixed (x,Q 2 ) 
values, there are drastic differences in the nucleon and 
in nuclei. For instance, resonance structure can be ob- 
served in the A resonance region in deuterium but not 
at all in iron. However, the effect of this redistribution is 
smaller in nuclei, such that the resonance region structure 
function ne arly reproduces the DIS structure function in 
nuclei [l4l Il8j |. and the ratios of the nuclear structure 
function in the resonance region rep roduce the observed 
EMC effect with high precision |4Cj . 

West points out the need to reconcile the difference 
between the fundamental asymptotic QCD sum rule, 

j\jF 2 A -±Fi>)dx = 0, (3) 

based on energy-momentum conservation, and the nomi- 
nal observation of the EMC effect that the nuclear struc- 
ture function is not simply A times that of a nucleon |ll| . 
The new data presented here (Fig.^J) indicate agreement 
with this sum rule already at the low Q 2 values here ob- 
served. As a quantitative example, the integral in Eq. |3| 
becomes of the order of 4 x 10~ 3 compared to individual 
moments of -0.2 at Q 2 = 2.9 GeV 2 . 

The moments of the structure function Fo can be de- 
termined theoretically on the lattice 0, 142ll43ll . While 
contributions from disconnected diagrams [4l| make it 
more difficult to calculate the separate proton and neu- 
tron moments on the lattice, these contributions cancel 
in the non-singlet combination M n (p) — M n (n). To com- 
pare our results with lattice calculations we extracted the 
difference between the proton and neutron moments for 
n=2 and n=4. We assume that the deuteron moment 
is equal to the sum of proton and neutron, and then 
determining the p — n moment from the proton |35j and 
deuteron moments, taking M n (p — n) = 2M n (p) — M n (d) . 
Because the proton and deuteron moments are sometimes 





This work 
Q 2 « 4 GeV 2 


Detmold 
et al. [421 


Dolgov 
et al. [411 


Gockeler 
et al. [43] 


n = 2 
n = 4 


0.049(17) 
0.015(03) 


0.059(8) 
0.008(3) 


0.269 
0.078 


0.245 
0.059 



TABLE III: Moments of the F2 structure function for the 
difference p — n. Experimental results for Q 2 « 4 GeV 2 from 
the present work are compared to lattice calculations at 4 
GeV 2 . 



extracted at slightly different Q 2 values, we combine our 
extracted deuteron moments with the nearest proton mo- 
ments, scaling the proton to the correct Q 2 value using 
the Q 2 dependence of the simple fit shown in Fig. [21 
The extracted values for the M2 and M4 moments for 
p — n are independent of Q 2 above 2 GeV 2 . The ex- 
perimental results shown in Table ITTT1 for Q 2 m 4 GeV 2 
come from combining the extracted values at Q 2 = 3.2 
and 4.5 GeV 2 , and are compared to lattice calculation at 
Q 2 = 4 GeV 2 . Because the proton and neutron n = 2 
moments are comparable in size, there is a large cancel- 
lation in the difference which leads to the large relative 
uncertainty. 

The n = 2 moment from Detmold et al. jl^ is m excel- 
lent agreement with the measured data. For n = 4, the 
small discrepancy between the lattice calculation and our 
experimental result could be due to higher twist effects, 
which are not included in the lattice result, although the 
Q 2 dependence of the moments does not indicate that 
these are large. In addition, no nuclear effects were taken 
into consideration when extracting the neutron moment 
from deuterium data. These effects seem to be small 
when averaged over the entire x range but they might 
still have some non-negligible contribution. It should also 
be noted that there are still open issues for lattice calcu- 
lations, such as chiral extrapolation, volume dependence, 
or renormalization. To demonstrate this, we also show 
the results of Dolgov et al. , and Gockeler et al. ^3 • 
The main difference between the lattice calculations pre- 
sented here is the chiral extrapolations used. In Ref. |41| . 
the lattice results are ex trap olated linearly to the phys- 
ical limit, while in Ref. [43, the extrapolation includes 
the correct chiral behavior from chiral effective theory. 

We note that comparisons between lattice and nomi- 
nal data formed from pdf-based fits have been performed 
previously [42j. We stress that such fits do not ade- 
quately account for the large x regime where they are 
unconstrained by data. Moreover, substantial uncertain- 
ties exist in the down quark distribution, d(x), associated 
with assumptions utilized in extracting neutron results 
from deuteron data, as well as the unknown behavior of 
d/u as x —* 1. 
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IV. CONCLUSIONS 

In conclusion, we utilized inclusive electron-nucleus 
scattering data to obtain nuclear structure function mo- 
ments for iron and deuterium. The new data are par- 
ticularly important for moment calculations at low Q 2 , 
where there was a paucity of previous data. Moreover, at 
low Q 2 and higher n, the need for large x data increases 
as this regime comes to dominate the moments. 

Negligible Q 2 dependence is observed in the lower order 
moments, indicating agreement with asymptotic predic- 
tions and minimal higher twist effects. This is surprising, 
given that the data extend to quite low Q 2 values. 

The n = 2 moment, related to the momentum sum 
rule, is here presented for both iron and deuterium. Ad- 
ditionally, a neutron momentum was formed by subtract- 
ing existing proton data from the deuterium data. The 
measured iron moments were found to agree with mo- 



ments simply constructed from these neutrons and pro- 
tons. This observation has interesting implications for 
interpretations of the EMC effect. 

Finally, these neutron and proton moment data allow 
for comparison with lattice QCD calculations. The ex- 
tracted non-singlet moments provide the first direct com- 
parison to lattice calculations of the non-singlet moments 
and the results are in good agreement with the calcula- 
tion of Ref. 0. 
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